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1. Introduction 



The theory of quarkonia production has recently undergone a number of significant 
developments. Progress in this area was stimulated by the discovery of Braaten and Yuan 
in 1993 that parton fragmentation represents the dominant source of prompt quarkonia at 
high energies These authors reahzed that the fragmentation functions which specify 
the probabihty for gluons and heavy quarks to hadronize into quarkonia bound states can 
be calculated starting from first principles in QCD. Over the past few years, the most 
important fragmentation functions for S and P-wave quarkonia have been computed to 
lowest order [^j-Q. They have subsequently been utilized to investigate a range of issues 
related to quarkonium phenomenology 

One of the most interesting applications of this fragmentation work has been to the 
study of prompt charmonium production at the Tevatron Comparison between 

theoretical predictions and experimental measurements of J/t/j and t/j' differential cross 
sections clearly demonstrate that fragmentation dominates over all other charmonium cre- 
ation mechanisms at large transverse momenta. Orders of magnitude discrepancies be- 
tween cross section computations and observations are significantly reduced when gluon 
and charm fragmentation processes are taken into account. The predicted prompt J/i(; 
rate then qualitatively agrees with recent CDF data [0. In the case of "0' production, the 
fragmentation results substantially improve upon earlier differential cross section predic- 
tions. But they still underestimate the number of '(/'"s observed at the Tevatron by more 
than a factor of 30. This large discrepancy between theory and experiment indicates that 
some important production mechanism beyond the simplest g ^ ij)' and c — > V' fragmen- 
tation processes needs to be included. A number of possible resolutions to the -0' problem 
have been suggested [p!0H T3|, and experimental studies which may confirm or reject these 
proposals are underway. 

Although the fragmentation picture of quarkonium production has provided valuable 
insight, its validity domain is restricted to high energies. The approximations that en- 
ter into fragmentation function computations break down when a quarkonium's energy 
becomes comparable to its mass. Fragmentation predictions for charmonium differential 
cross sections are therefore unreliable at low transverse momenta. Similarly, fragmentation 
results for bottomonium production are untrustworthy throughout the p± < 15 GeV re- 
gion where Tevatron data exists. Yet even in this low energy regime, Upsilon cross section 
measurements significantly disagree with existing theoretical predictions . 
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The discrepancies in the Upsilon sector provide the primary motivation for the work 
which we present in this paper. Building upon several recently developed ideas in heavy 
quarkonium physics, we identify a large class of color-octet diagrams that mediate quarko- 
nia production at all energies and reduce to the dominant set of gluon fragmentation graphs 
in the high p± limit. These diagrams enter at the same order into direct production pro- 
cesses as previously considered color-singlet graphs and must be included for consistency. 
As we shall see, the addition of color-octet contributions to Upsilon production computa- 
tions significantly improves agreement between theory and experiment. 

Our article is organized as follows. In section 2, we discuss the short and long distance 
factorization of quarkonia production amplitudes. We then calculate the leading color-octet 
contributions to quarkonia differential cross sections. In section 3, we apply the color- 
octet cross sections to the study of Upsilon and Psi production at the Tevatron. Using 
CDF data, we extract numerical values for bottomonia and charmonia color-octet matrix 
elements which we compare with NRQCD scaling relations in section 4. We next analyze 
the spin alignment of quarkonia generated via the color-octet mechanism as a function of 
transverse momentum in section 5. Finally, we summarize our findings in section 6 and 
close with some thoughts on directions for future work. 



2. Color octet contributions to quarkonia cross sections 

The physics of quarkonia involves several different energy scales that are separated by 
the small velocity v of the two heavy constituents inside QQ bound states [jl5| . The first is 



set by the quark mass Mq which fixes the distance range for QQ creation and annihilation 
processes. The second is given by the heavy constituents' typical momentum MqV which is 
inversely proportional to the bound state's spatial size. The kinetic energy Mqv"^ represents 
a third important quarkonium energy scale and determines the constituents' interaction 
time. As Mq — ^ oo, the heavy constituents' velocity v ~ 1/ log Mq tends towards zero, 
and these three scales become widely separated: 

{MQV^f^{MQvf<^Ml. (2.1) 

This hierarchy is well realized in the real world for bottomonia with v'^ ^ 0.08 and to a 
lesser extent for charmonia with v'^ ~ 0.23. 
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The presence of these different distance scales renders the study of quarkonia both in- 
teresting and chaUenging. An effective field theory formalism called Nonrelativistic Quan- 
tum Chromodynamics (NRQCD) has been established to systematically keep track of this 
scale hierarchy |1^,16[. It is based upon a double power series expansion in the strong 



interaction fine structure constant ag = g1/^T^ and the small velocity parameter v. Short 
distance processes which occur at energies comparable to or greater than the heavy quark 
mass are perturbatively computable within the effective theory. Long distance effects are 
consolidated into nonperturbative matrix elements whose values must be determined from 
experiment or the lattice. High and low energy interactions are thus cleanly separated in 
NRQCD, and simple power counting rules methodically determine the dominant contribu- 
tions to various quarkonium processes [0. We will work within the NRQCD framework 



throughout this paper. 

Before the development of NRQCD, quarkonia were generally treated as nonrelativis- 
tic bound states of a heavy quark and antiquark in a static gluon field that sets up an 
instantaneous confining potential [jl^ . Although this picture has enjoyed a remarkable de- 
gree of phenomenological success, it fails to take into account gluons inside a quarkonium 
with wavelengths much greater than the bound state's characteristic size. The presence 
of such low energy gluons implies that the heavy quark and antiquark cannot always be 
regarded as residing in a color-singlet configuration. 

This shortcoming of the potential model approach is rectified in NRQCD. Dynamical 
gluons enter into Fock state decompositions of physical quarkonium states. For example, 
wavefunctions for S-wave orthoquarkonia schematically appear in Coulomb gauge as 

I = I QQ ) + 0{v) I QQ ['Pf] g) 

+ 0{v') I QQ fs['''^] gg) + 0{v') \ QQ [^^J^^] g) + 0{v') \ QQ fD^J^'^] gg) + ■ ■ ■ . 

(2.2) 

The angular momentum quantum numbers of the QQ pairs within the various Fock com- 
ponents are indicated in spectroscopic notation inside the square brackets, and their color 
configurations are labeled by singlet or octet superscripts. S-wave quarkonia reactions 
proceed at lowest order in the velocity expansion through the first Fock state in ( |2.2| ). The 
QQ pair in this 0{1) Fock component has precisely the same angular momentum and color 
quantum numbers as the full quarkonium state. The NRQCD description of i/jq production 
or annihilation thus reduces to the familiar color-singlet model in the v limit. 



P-wave orthoquarkonia processes are more complicated. They are mediated at leading 
order in the velocity expansion by the first two Fock states in the Xqj wavefunction 

I x«.) = I QQ ) + 0{v) I QQ [^Sf^] g) + .-.. (2.3) 

Contributions to Xqj formation which involve the first color-singlet Fock component are 
proportional to the squared derivative of the P-wave bound state's wavefunction at the 
origin. This quantity counts as 0{v^) in the NRQCD velocity expansion. On the other 
hand, short distance creation of the S-wave color-octet QQ pair in the second Fock state 
of (|2.3|) takes place at 0{v^). The subsequent long distance evolution of the QQ^S'f'\ 



pair into a colorless Xqj hadron via the emission or absorption of a soft gluon costs an 
additional power of v in the amplitude and v"^ in the rate. So Xqj production proceeds at 
O(f^) in both the color-singlet and color-octet mechanisms [Q. 

It is useful to draw a picture which clarifies the difference between these two quarko- 
nium production mechanisms. In fig. |^a, we illustrate a typical color-singlet Feynman 
graph that creates a Xqj- The 0(ag) hard scattering forms a colorless QQ'^P^^] pair at a 
short distance scale. The heavy quark and antiquark fiy out from the initial collision point 
in nearly parallel directions and almost on mass shell. After exchanging many soft gluons, 
the color-singlet QQ pair eventually hadronizes at a long distance into a Xq.i bound state. 
In fig. |l]b, the 0{a^g) high energy collision creates a QQ pair in a color-octet configuration. 
Far away from the collision point, the QQf'^S'^^''] pair emits a long wavelength gluon which 
bleeds off its color but carries away virtually no energy or momentum. The heavy pair 
thus transforms into a colorless Xqj quarkonium state. 

Color-singlet quarkonium production has been studied for years. Leading order color- 
singlet differential cross sections were calculated a decade ago [^9|-^, and some total 



cross section formulae have been evaluated at next-to-leading order as well In con- 
trast, color-octet contributions to quarkonium processes have only recently begun to be 
considered. While the latter mechanism may be less familiar than the former, it is cer- 
tainly not less important. As can be seen in fig. |l|, color-singlet and color-octet graphs arise 
at the same order in a^. Moreover, the color-octet diagrams dominate at high energies. 
So a complete analysis of quarkonium production at hadron colliders must include both 
mechanisms. 
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Color-octet creation of S-wave orthoquarkonia starts at 0{al) via the qq i/jq Feyn- 
man diagram shown in fig. As indicated by the dashed hne in the figure, the total 
color-octet amplitude factorizes into short and long distance pieces: 



A{qq ^ MnS))^^^^^ = A{qq ^ QQp^r^])short distanceX^(QQp^r^] ^ ^Q(n5))i„„gdistanc 

(2.4) 

The short distance part of the quark-antiquark fusion process lies to the left of the dashed 
line in fig. ^ and can be calculated within perturbative QCD. On the other hand, the long 
distance component to the right of the dashed line involves nonperturbative physics. We 
will consider these two factors in turn. 

To begin, we examine the subamplitude for virtual gluon formation of a color-octet 
QQ pair 

AqI ^ QQ?sf\)^ = E Y^H^+q-.si)gsl^{Ta)]v{^-q;s2){ls,; ^.2|1^.) (3i; 3i|86) 

(2.5) 

which is illustrated in fig. ^. The heavy quark and antiquark propagate nearly on-shell 
with a large combined momentum P and a small relative momentum q. The spin and color 
quantum numbers for the QQ pair are projected out by the sums over the SU (2) and SU (3) 
Clebsch-Gordan coefficients {^si; ^S2\lSz) and (3z;3j|86) = -\/2(Tb);^. The subamplitude 
may be rewritten in the trace form 

A{g: ^ QQ['si%)^ = V2g,Tr{TMTr[j^P,s. (P; q)] (2.6) 

where 

Pc^c: (P:n) = \^ v( . 

2 - V 2 ^^'2 ' 2 

Sl,S2 

denotes a spin projection operator. This 4x4 matrix reduces to the covariant expressions 

(2 8) 

Pis. {P; q) = (| - ^ - Mc,)r (P; s.) (| + ^ + M^) 



Pss. (P; q)=J2v{^~q; S2)u{^ + q; Si) (^Si; ^Ssl^^.) (2.7) 



^ Spin projection matrices were originally introduced in refs. |23,^|. The u and v spinors in 
our definition of Pss. correspond to an outgoing QQ pair. They differ from the v and u spinors 
appearing in these earlier articles' spin projection operators. Our (+, — , — , — ) metric signature 



convention is also opposite to that of ref. |24]. 
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up to 0{q'^) corrections. The projection operators' dependence on q is irrelevant for S-wave 
QQ production. The g* —>■ QQ[^S[^'^] subamphtude consequently simplifies to 

Aigl - QQ['S[%)^ = gsMe;{P; S,)SaB (2.9) 

where M = 2Mq denotes the mass of the heavy quark-antiquark pair. Combining this 
result with the remaining perturbative part of the Feynman graph in fig. we determine 
the short distance matrix element in eqn. (p.4| ): 

2 

A{q{pi)q{p2) ^ QQ['^r^]a)short distance = J^v{P2h''TaU{pi)e*^{pi + P2; S . (2.10) 

Once the color-octet QQ pair has been formed, it can evolve into a physical i/^q bound 
state in many different ways. For example, it may emit two long wavelength gluons in 
a double chromoelectric dipole transition and turn into a color-singlet QQ pair that has 
order unity overlap with the full i/jq wavefunction in ( |2.2| ). Alternatively, it may absorb 
a soft gluon and transform into a Xqj through the 0{v) Fock state in (|2.3|) . The P-wave 
quarkonium can subsequently emit a photon and decay to the S-wave state. All possible 
ways to produce a i/jq starting from a QQ[^s[^^] pair are encoded into the long distance 
amplitude in eqn. (|2.4| ). This factor is pictorially represented by the big blob in fig. ^ 
with the long wavelength gluon tail. It cannot be computed from first principles without 
resorting to lattice QCD. But since the low energy dependence of interactions that take 
place inside the blob is negligible, the nonperturbative amplitude is simply a constant: 

A{QWS['^] - V'«(n5)) longdistance = ^Onst. (2.11) 

The total color-octet contribution to V'q production thus involves one free parameter which 
we will eventually have to fit. 

We turn now to investigate the gluon fusion process gg i^q- The color-octet ampli- 
tude for this reaction again factorizes into short and long distance pieces: 

A{gg ^ V'Q(^^))„,tet = ^(99 - QQP^f ^])shortdistanceX^(QQp^f ^ Mu S)) ^^^^ ^.^^^^ 

(2.12) 

The long distance term is precisely the same as in eqn. ( p.ll|) . The short distance factor 
involves at lowest order the three Feynman diagrams displayed in fig. ^. Evaluating the 
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high energy gluonic matrix element in a similar fashion to its quark counterpart in ( |2.10|) , 
we deduce 

• 2 

A{ga{pi)9h{p2) ^ QQ['^f^]c)short distance = fahcS^ {pi)e'' {p2)e%Pl + S 

2 2 (2.13) 

^ 2 / 2 ^^2 [(^^2 - Pl)ag^lu + 'i{Plug^la " P2m^!^<t)] • 
Pi + P2 -''^ 

It is important to observe that this expression vanishes when both incoming gluons are on- 
shell. This result is consistent with Yang's theorem which forbids a massive J = 1 vector 
boson from decaying to two massless J = 1 bosons [^. In fact, the theorem requires 
that the on-shell gg — > iI)q amplitude vanish to all orders. The color-octet diagrams in 
fig. ^ consequently participate in quarkonium production only as subgraphs inside more 
complicated reactions. 

0{a1) scattering processes create quarkonia at small transverse momenta that differ 
from zero only as a result of the incident partons' intrinsic hadronic motion. Since outgoing 
states from 2^1 collisions are lost down the beampipe, we need to consider 2 2 tran- 
sitions which produce quarkonia with experimentally measurable p±. Such reactions start 
at 0{al) and proceed in both the color-singlet and color-octet mechanisms through the 
partonic channels qq QQ[^^+^Lf^]g, gq QQ[^^+^Lf^]q and gg QQ[^'^+^Lf^]g. 



The differential cross sections for the color-singlet modes were calculated in refs. [[IQl -pl 
We present their color-octet analogues below. 

The Feynman diagrams which generate i/'q quarkonia in the three different partonic 
modes through the color-octet mechanism are illustrated in figs. 5a, 5b and 5c. The shaded 
circles in these figures represent the qq — > if^Q and gg if^Q amplitudes in eqns. ( p^ ) and 
( p.l2| ). The sums of the graphs in each channel yield gauge invariant amplitudes. Their 
spin and color averaged squares are given by 

E \A{qq Ml' = 27M2 iu{s-M^r ^ ^ ^ ^ " H-^sIV'q) (2.14a) 

E \A{gq - = - g^,- ga(f_M2)2 [^^^ + ^ ) - H-^8(^g) (2.146) 

^ 2_ 327r3a3 27 {si + iu + us) - IdM"^ 

X^\A[gg-.i;,g}\ -- ^ M^)it - M^)iu - M^f 

X [{P + iu + u^f - M^{i + u){2p + iu + 2v?) + M^{P + iu + u^y^Msii^q) (2.14c) 

where Ais{ipQ) = \A{_QQ[^s[^'^] — > 'i/'Q(n5')) p. The squared matrix elements in ( |2.14a| ) 
and (|2.146| ) were determined using standard spinor summation techniques. The result in 
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( p.l4cj) was obtained from a simple helicity amplitude method which obviated the need 
to compute a large number of interference terms. We calculated the 24 different helicity 
amplitudes for gg —>■ t/^gg scattering by choosing explicit representations for the mass- 
less and massive vector boson momenta and polarizations in the parton center-of-mass 
frame, i The amplitudes were then individually squared and added together using the 
high energy physics Mathematica package FEYNCALC [ |26[ |. 

The squared matrix element in eqn. ( p.l4a| ) diverges when either of the partonic 



Mandelstam invariants t or u approach zero. The source of these long distance infinities can 
be identified by looking at fig. ||a. As the outgoing gluon in a q'^ ^ i/^gg collision becomes 
collinear with either the incoming quark or antiquark, the exchanged fermion in the t or 
M-channel graphs goes on shell and its propagator blows up. The divergence generated by 
the collinear splitting of the incident parton may be factored into its distribution function. 
The hard scattering reaction then reduces to the qq tpQ Born level process illustrated in 
fig. ^. A similar w-channel collinear divergence develops in the squared gq tpqq amplitude 
when the four-momenta of the incoming gluon and outgoing quark become proportional. 
In contrast, the t-channel diagram in fig. jp) remains finite. The pole in the intermediate 
propagator is canceled by the zero in the gg if^Q amplitude as the t-channel gluon goes 
on shell. This same cancellation renders ^ \ A{gg ?/'q(7)P free of long distance infinities. 

The squared color-octet amplitudes in ( |2.14| ) enter into the partonic differential cross 
section 

-^(a6 V'Qc)octet = '-^(^^ ^ V'qc) P (2.15) 

which appears in turn within the triple differential hadronic expression 

, ^ ^ , {AB ^ iIJqX) octet = '2p±y^^XaXbfa/A{Xa)fb/B{Xb)^{ab ^ ijQc)octet- (2.16) 

dy3dy4:dp± ^ dt 

aoc 

The partonic cross section is multiplied by distribution functions fa/Ai^a) and fb/sixb) 
that specify the probability of finding partons a and b inside hadrons A and B carrying 
momentum fractions Xa and Xb- The product is then summed over parton channels. The 
resulting hadronic cross section is a function of the ifjg and recoiling jet rapidities and 
1/4 and their common transverse momentum p± . 



^ Parity and crossing symmetry relations imply that only 4 of the 24 gg ipgg helicity 
amplitudes are independent. 



8 



It is instructive to examine the high energy hmit of color-octet quarkonia production. 
As the partonic Mandelstam invariants grow to infinity, the cross section in eqn. ( |2.15| ) 
reduces to 



^(a6-.^^c)octet ^{ab^g*c)x^^y x\A{g*^i;c,)\'. (2.17) 



This asymptotic expression has a simple gluon fragmentation interpretation. The first 
factor represents the differential cross section for producing a high energy virtual gluon. 
The second term comes from the square of the gluon's propagator. The last factor equals 



the square of the amplitude in eqn. (|2.9| ) times Aisi^JQ) and determines the virtual gluon's 
probability to hadronize into a i/'q bound state. The gluon fragmentation picture for heavy 
quarkonium production is thus precisely recovered in the high energy limit |jl[]. 

Gluon fragmentation via the color-octet mechanism represents the dominant source of 
large p± quarkonia at hadron colliders P,|6|,p!0[] . The total cross section for i/jq production 
reduces at high energies to the fragmentation form 



dysdy^dp 



± 



(2.18) 

The gluon fragmentation function evaluated at the factorization scale fi = M is readily 
identified from eqn. ( p.l7| ): 

Z>,,,„(..M) = ^-°-<^;^-'^^) .(l-.). (2.19) 

Leading log QCD corrections to this result may be summed up using the Altarelli-Parisi 
equation 



dDg^i,Q as[ix) f dy ,z . 



where 



[2.21) 



denotes the gluon splitting function for n/ active quark fiavors. At high energies, the 
fragmentation approximation in ( |2.18| ) incorporates sizable 0(log(i?^/M^)) renormaliza- 
tion effects, and its intrinsic 0{M'^ / E'^) errors are negligible. In contrast, the color-octet 
formula in (|2.16| ) does not include any QCD corrections which are small at low pj_, but it 
retains full dependence upon all 0{M'^ /E'^) terms. These two forms for the i/jq differential 
cross section are thus complementary. 
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We may exercise our perturbative freedom to supplement the 0{a^) color-octet cross 
section with the leading logarithms from the gluon fragmentation approach. We implement 
this choice as follows: 



\ T w J rrag without QOD running 

(2.22) 

The ratio on the RHS of this hybrid expression reduces to unity at low energy, whereas the 
first factor over the denominator approaches unity at high energy. Eqn. (|2.22| ) therefore 
smoothly interpolates between the two asymptotic limits in eqns. ( |2.16| ) and ( |2.18| ). We 
will use this final color-octet formula to study bottomonia and charmonia production at 
the Tevatron in the following section. 



3. Upsilon and Psi production at the Tevatron 

The CDF collaboration has recently measured Upsilon production at the Tevatron 



for the first time Their reported T(15'), T(25') and T{3S) transverse momentum 

differential cross sections seriously disagree with predictions based upon 0{al) color-singlet 
cross section formulae. While some discrepancy between the theoretical and experimental 
findings may be attributed to various parameter uncertainties and higher order corrections, 
the magnitude of the disparity strongly suggests that additional production mechanisms 
beyond the simplest color-singlet processes are at work. As we shall see, inclusion of 
color-octet reactions brings the theoretical cross sections into line with the Upsilon data. 

We first investigate the differential cross section for T(15') production in the rapidity 
range \y\ < 0.4. The color-singlet distribution is illustrated by the dashed curve in 
fig. la. I It includes direct T(15') production as well as radiative feeddown from Xbj{^P) 
and XbjC^P) states. Additional contributions from strong and electromagnetic decays of 
T(25') and T(35') are negligible. In computing the color-singlet cross sections, we adopted 
the radial wavefunction values tabulated in ref. for the Buchmiiller-Tye QCD potential 
||2^ and the associated bottom quark mass Mi, = 4.88 GeV. The color-singlet prediction 
for T(IS') production clearly underestimates the CDF data points shown in the figure. 

The dot-dashed curve in fig. |^ displays the color-octet contribution to the T(IS') 
differential cross section. The best fit value M8{T{1S)) = (3.66 ± 0.27) x lO'^ GeV^ 



^ The differential cross sections displayed in fig. ^ and subsequent figures were calculated using 
the MRSDO parton distribution functions. 
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for its squared long distance matrix element determines the probability for a 66['^5'{^''] 
pair to hadronize into an T(15'). Following eqn. ( p.22| ), we have incorporated leading 
log QCD corrections into the color-octet cross section. The dot-dashed curve therefore 
approaches by construction the dotted gluon fragmentation curve based upon the same 
value for A^8(T(15')) at high p±. The two cross sections diverge, however, at transverse 
momenta comparable to the T(IS') mass. This disparity demonstrates the breakdown of 
the fragmentation picture in the low energy regime. 

The sum of the color-singlet and color-octet transverse momentum distributions is 
illustrated by the solid curve in fig. ^ It represents the total theoretical T(15') differential 
cross section. The solid curve's shape fits the data quite well except at the lowest point. 
At very small p_L, the color-singlet and color-octet cross sections are corrupted by coUinear 
divergences as discussed in section 2. Since the intrinsic motion of incident partons inside 
colliding hadrons renders the differential cross section uncertain for p±<2 GeV, we have 
not attempted to remove the divergences from the dashed and dot-dashed curves. Instead, 
we simply regard the portion of the plot which runs below p± = 2 GeV as untrustworthy. 
Above this low momentum cutoff, the predicted T(15') distribution matches the data. 

Analogous theoretical and experimental differential cross sections for T(2S') produc- 
tion in the rapidity region \y\ < 0.4 are displayed in fig. |^. The dashed color-singlet curve in 
the figure incorporates direct T(2S') production and radiative feeddown from Xb.i{2P). The 
dot-dashed color-octet and dotted gluon fragmentation curves correspond to the squared 
matrix element M8{T{2S)) = (2.25±0.38) x lO'^ GeV^. The sohd curve again illustrates 
the combined color-singlet and color-octet contributions to the T(2S) differential cross 
section. It clearly agrees with the CDF measurements much better than the color-singlet 
distribution alone. 

We consider next the T(3S') differential cross section for which data is limited but 
interesting. The extent to which radiative Xbj(3P) T(35')7 transitions mediate T(35') 
formation is unknown, for the n = 3 P-wave states have never been observed. Although the 
T(45') can electromagnetically decay to them, it lies above the BB threshold. As a result, 
Br(T(4S') X6j(3P)7) is very small, @ and the n = 3 P-wave wave levels are rarely pop- 
ulated at e+e~ colliders. But at the Tevatron, x&j(3-P) bottomonia are directly created in 
hadronic reactions. Since the integrated T(3S') cross section measured by CDF is 40 times 
larger than the color-singlet model's prediction when only T(35') production is taken into 



^ CLEO has recently set a 5% upper limit on non- BB decays of the T(45) 
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account, it is believed that n = 3 P-wave levels reside below the BB threshold and decay 
at a significant rate to the n = 3 S-wave state WM- We consequently include their contri- 
butions into the dashed color-singlet curve in fig. |^ using the measured n = 2 branching 
ratios Br(xbj(2-P) — > T(2S')7) as estimates for their unknown n = 3 counterparts. The dot- 
dashed, dotted and solid curves in the figure represent the color-octet, gluon fragmentation 
and total differential cross sections based upon 7Vl8(T(3S')) = (0.53 ±0.24) x 10""^ GeV^. 
Given all the uncertainties, we can only say that agreement between the few measured 
T{3S) data points and the theoretical prediction appears adequate. 

We now turn to the charmonium sector and examine prompt J/i(; production at the 
Tevatron within the pseudorapidity interval |r7| < 0.6. The J/tl^ cross section multiplied by 
the muon branching fraction Br (J/?/' = 0.0597 is illustrated in fig. ^. The dashed 

curve in the figure depicts the color-singlet differential cross section that includes direct 
J/ifj production and radiative feeddown from Xcj{^P) states. In calculating this curve, we 
employed the Buchmiiller-Tye charmonium wave function values at the origin tabulated 
in ref. and the associated charm quark mass = 1.48 GeV. Additional color-singlet 
contributions to the J/t/j cross section from i/;' production are small. The dot-dashed 
curve in fig. ^ shows the color-octet J/t/j cross section for the squared matrix element value 
MsiJ/il^) = (0.68 ± 0.03) X 10-3 GeV^. It is almost identical to the sohd curve which 
represents the sum of the color-singlet and color-octet differential distributions. Comparing 
the dot-dashed color-octet and dotted gluon fragmentation curves with the CDF data 
in fig. ^ we see that the fragmentation approximation appears to fit the experimental 
measurements better at small transverse momenta. But we must stress that this low 
energy agreement is fortuitous. The fragmentation prediction looks better only in the 
region where it becomes unreliable. 

We finally consider prompt t/^' production at the Tevatron. The observed differential 
cross section is 30 times larger than recent theoretical predictions based upon c ^ ij)' and 
g —* i/j' fragmentation calculations These two 0{v^) fragmentation processes create 

color-singlet cc['^S'^^'*] pairs at short distance scales and yield comparable i/j' production 
rates. At higher order in the NRQCD velocity expansion but lower order in perturbative 
QCD, the color-octet mechanism also mediates g —>■ fragmentation through an interme- 
diate cc['^S'^^'*] pair. The additional color-octet contribution to the i/j' cross section which 
was not previously taken into account can readily resolve the factor of 30 discrepancy 



between theory and data [12,13 
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We plot in fig. ^ the ij)' differential cross section multiplied by the muon branching 
fraction Br('i/'' — fx~^n~) = 0.0077. A pseudorapidity cut < 0.6 has been applied to the 
distributions shown in the figure. The dashed curve represents just direct t/j' production, 
for no other charmonium level below the DD threshold can decay to the n = 2 state. This 
color-singlet prediction drastically underestimates the measured differential cross section. 
Inclusion of color-octet contributions with A^g^i/j') = (0.20 ±0.01) x 10~^ GeV^ brings the 
magnitude and shape of the theoretical distribution into agreement with the data points 
displayed in fig. |TD|. The color-octet mechanism thus appears to solve the CDF t/j' problem. 



4. NRQCD matrix elements 

We have seen that heavy quarkonia production involves short and long distance 
physics. The high energy creation of color-singlet and color-octet QQ pairs is pertur- 
batively computable within the NRQCD effective theory. The subsequent low energy 
hadronization of these pairs into physical bound states is described in terms of nonrenor- 
malizable operator matrix elements. NRQCD scaling rules determine the relative impor- 
tance of different long distance matrix elements and yield relations among them. We will 
use these rules to check the consistency of the fitted Upsilon and Psi amplitude values 
which we obtained in the preceding section. 

We first recall the relations 

(0|O^e^0|0) = ^|i?(0)P = O(M^^^) 
{0\O^^\'Pj)\0) = -^{2J+1)\R'{0)\' = 0{Mlvl) 

between quarkonia radial wavefunctions at the origin and matrix elements of certain color- 
singlet four-quark operators in the NRQCD Lagrangian [|T^ . The details of these operators' 



definitions are discussed in ref. [|T5[, but they are not important to us. Instead, we are only 
interested in their scaling dependence upon the heavy quark mass Mq and velocity Vq. 
Using the Buchmiiller-Tye wavefunction information tabulated in ref. ||2^ , we list the 



numerical values of color-singlet matrix elements that are relevant for Upsilon and Psi 
production in Table I. 
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dnlnT-Sinf?lpt 


TvJnmPT'iral 


NROCD 


Matrix Element 


Value 


Scaling Order 


\ '-'1)1^/ 


3 9 X 10~^ GeV"^ 




(OlO^^H^p )|Q\ 


3.2 X 10"^ GeV^ 


c c 


(Olof S^)\0) 


2.5 X 10"^ GeV^ 


c c 




3.1 GeV^ 


-'"6 ^6 


/Q|(-)Xbi(iJ')/3D ^|Q\ 


6.1 GeV^ 


-'"6 




1.5 GeV^ 


Muvl 




7.1 GeV^ 






1.2 GeV^ 






7.7 GeV^ 





Table I 

The entries in this table are consistent with the NRQCD power counting rules. In the 
— > limit, the long distance quarkonia matrix elements would be independent of the 
radial quantum number n. In actuality, the n dependence of the hb sector entries is 
generally smaller than that of the cc values. This trend simply indicates that the NRQCD 
velocity expansion works better for bottomonia with ~ 0.08 than for charmonia with 
vl ~ 0.23. 

NRQCD scaling rules can also be applied to color-octet four quark operators. We 
focus upon the matrix elements 

{Q\0t''{'S,m = 0{Mlvl) 



{^\Ol''\^S,m = 0{Mlv%) 



(4.2) 



14 



that determine the probabihty for QQ[^s[^^] pairs to hadronize into tpQ and Xqj quarkonia. 
These matrix elements appear in the decomposition of the total long distance squared 
amplitude 

MsiMn)) = ^ 5^ {(OlO^'^^-^^e^OlO) Bv{Mm) ^ ^^(n)) 

J=0,l,2 

(4.3) 

which reflects the different ways for a color-octet pair to materialize as a final i(jq. The 
infinite number of terms represented by the ellipsis contribute at higher order in the Vq 
or as expansions to color-octet i/'q production than those explicitly displayed. So we will 
neglect them for simplicity. 

We can extract numerical values for the NRQCD color-octet matrix elements in 
eqn. ( ^.2p by inserting our fitted values for M.s{i/^q) into eqn. ( |4.3D . Combining M.^[J /^) = 
6.8 X 10"^ GeV^ and MsW) = 2.0 x lO'^ GeV^ with the recent CDF finding that 32.3% 
prompt J/i/^'s originate from Xcj decays ||3^, we deduce the charmonium color-octet matrix 



elements 

{0\O^^'^CSi)\0) = 1.2 X 10"^ GeV^ (4.4a) 

{0\of{^Si)\0) = 7.3 X 10-3 GeV^ (4.46) 

{0\O^^'CSi)\0) = 1.6 X 10-2 GeV^ (4.4c) 

In the absence of additional bottomonium experimental information, we employ the 
NRQCD relation 

/0|oJ(-5)^3, 



(4.5) 



to scale the Psi values in ( [4.4q| ) and ( [1.45|) to the Upsilon sector. We then obtain the 



color-octet matrix elements listed in Table II. 
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Color- Octet 


TMiimeriral 


NROCD 


Matrix Element 


Value 


Scaling Order 




1 2 X in~^ GeV"^ 






1.6 X 10"2 GeV^ 


c 


(0\of Si)\0) 

\^\^8 \ 


7 3 X 10"^ GeV^ 






9.5 X 10"^ GeV^ 


M?vl 




4.3 X 10-1 GeV^ 


MM 




5.7 X 10-3 GeV^ 


MM 




5.2 X 10-1 GeV^ 


MM 



Table II 

The numerical entries in this second table should be regarded as approximate esti- 
mates. Nonetheless, we see that the NRQCD scaling rules hold reasonably well for the 
color-octet matrix elements. The magnitudes of the entries in Table II appear consistent 
with each other and with those in Table I. These results thus provide a useful check on 
the entire color-octet quarkonia production picture. 



5. Quarkonia spin alignment 

We have so far not considered final state spin information in our analysis. Previous 
investigations have found that quarkonia generated at large transverse momenta via gluon 
fragmentation processes like gg gg* ggXqj cire significantly polarized [jl^. The source 
of this spin alignment can be traced to the fragmenting gluon. At high energies, the virtual 
g* is off-shell only by 0{M'^ / E"^) and is therefore nearly real and transverse. The <5Q['^5'j^''] 
pair which the gluon forms at short distance scale must inherit its transverse alignment 
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in order to conserve angular momentum. The subsequent long distance evolution of the 
color-octet pair into a physical quarkonium state preserves the spins of the heavy quark 
and antiquark up to 0(Aqcd/-^) corrections. The final quarkonium's polarization state 
is then fixed by heavy quark spin symmetry ||12|| . 

It is interesting to examine the spin alignments of quarkonia produced via the color- 
octet mechanism at low energies as well. To do so, we decompose the partonic squared 
amplitudes in eqn. (|2.14[) into separate longitudinal and transverse components: 



h=0 ^ ' 

I-4(m-QQ' - 27M2 (i-Af2)4 W + n'')^tu] (6.16) 



h=0 



\h\=l 



X [4(s2 + u^) - su] (b.ld) 
2_^\M99^QQ [S\']9)\ = :7T7^J^^J72Y2^^ + u )tu 



h=0 



27( si + iu + us) - 
[(s-M2)(£-M2)({i-M2)]' 



|h|=i ^ ^ 

27isi + iu + us) -19M^ 

^ 1 ( 5 1 / ) 

[(s-M2)(t-M2)(w-M2)]^ 

The results in the qq — > QQ''^'' ^ S^^]g and (/g — QQ^^'^ ^ s'f''^]q channels were calcu- 
lated using the covariant longitudinal and transverse spin sums discussed in ref. . The 
gg QQ^ S^^]g squared matrix elements were simply derived from our previously de- 
scribed helicity amplitude method. Summing over helicities in each mode, we reproduce 
the unpolarized expressions in eqn. ( |2.14| ). 

The longitudinal squared amplitudes in ( |5.1aD , ( |5.1cj ) and ( p.leD vanish as s/M"^ oo, 
whereas their transverse counterparts in (|5.16| ) , ( |5.1cj ) and (|5.1/| ) approach constant limits. 
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The gluon fragmentation prediction of total transverse alignment is therefore recovered at 



high energies. This asymptotic behavior can be seen in fig. 11 where we plot the ratio 



-'•'trans — i W--^) 



h 



of transverse and total 66["^S'^^''] and ccf'^S'^^''] differential cross sections as a function of 
p±. i Fig. 11 also demonstrates that color-octet pairs are strongly transverse even at 
low energies. This surprising phenomenon has a simple explanation. In the dominant 
99 ~^ QQ^'^\^s[^^]g channel, longitudinal color-octet pair production along the beam axis 
is prohibited by angular momentum conservation. The squared h = amplitude in ( p.l^ ) 
therefore vanishes as t = — — M^)(l — cos6'*) or u = — — M^)(l + cos6'*) approaches 
zero, while its \h\ = 1 counterpart in ( 5.1/1) remains finite. Since small angle scattering 



dominates this 2^2 partonic process, most QQ['^S[^^] pairs are formed at | cos 6**1 ^ 1. 
As a result, they are transversely aligned independent of their energy. The ratios of the 
h = and \h\ = 1 squared amplitudes in the qq QQ^'^\'^s[^'^]g and gq QQ^'^\'^s[^^]q 
channels similarly vanish as t — or w — > 0. 

The nearly 100% transverse alignment of cc[^5'(^'*] pairs should be directly observable 
in '0' production at the Tevatron if the color-octet mechanism resolves the i/;' problem. 
Color-octet production then dominates over all other prompt t/j' sources by more than a 
factor of 30. Heavy quark spin symmetry guarantees that '0' spin alignments equal those of 
their cc['^5'P'*] progenitors in the oo limit fl^. Finite charm mass effects somewhat 

degrade final "0' polarizations. But experimental observation of a large transverse '0' spin 
alignment would provide strong support for the color-octet production picture. 

Other S and P-wave charmonia and bottomonia can also inherit a sizable spin align- 
ment. Their helicity level populations depend upon the color-octet matrix elements in 
Table II. Detailed predictions for quarkonium polarizations resulting from color-singlet 
and color-octet production will be presented elsewhere |]31|. But we emphasize here that 
experimental measurements of Upsilon, Psi and Chi spin alignments will provide valuable 
tests of our understanding of basic quarkonium physics. 



^ The shapes of the 66['^>S'|^^] and cc[^S'^^^] curves in fig. 11 are not exactly the same since their 



respective \y\ < 0.4 rapidity and \r]\ < 0.6 pseudorapidity cuts are different. 
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6. Conclusion 

In this paper, we have examined color-octet production of heavy quarkonia at hadron 
colliders. Color-octet diagrams arise at the same order in perturbative QCD as their 
color-singlet analogues and reduce to the dominant set of gluon fragmentation graphs in 
the high energy limit. They must therefore be included in bottomonia and charmonia 
production computations. As we have seen, color-octet contributions to Upsilon and Psi 
differential cross sections eliminate large disparities between earlier predictions and recent 
measurements. The long distance matrix element values needed to bring theory into line 
with experiment are consistent with NRQCD scaling rules. We have also found that the 
color-octet mechanism yields QQ[^S[^^] pairs which are nearly 100% transversely aligned 
at all p_i. Observation of large quarkonia spin alignments would provide strong support 
for the color-octet production picture. 

A number of extensions of this work would be interesting to pursue. For example, a 
low transverse momentum study should reveal that quarkonia distributions vanish rather 
than diverge as pj_ 0. In order to compute the cross section turnover, it will be necessary 
to factorize coUinear singularities in both the color-singlet and color-octet channels into 
incident parton distribution functions. This factorization also needs to be performed before 
differential quarkonia distributions can be integrated. We expect color-octet contributions 
to significantly enhance existing predictions for total Upsilon and Psi production rates. 
Comparison of integrated quarkonium cross sections with fixed target measurements will 
provide an important means for testing the long distance matrix element values which we 
have found using collider data. Finally, all the results in this article can be applied to the 
study of quarkonia production at the next generation of hadron colliders. We look forward 
to further surprises in quarkonium physics coming from machines like the LHC well into 
the next century. 
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Figure Captions 



Fig. 1. Typical Feynman grapiis wliicli contribute to Xqj production at 0{a^g) through 
the (a) color-singlet and (b) color-octet mechanisms. The short-distance scatter- 
ing collisions respectively create QQ pairs in ^P^^ and ^S^^ configurations. The 
QQ pairs hadronize at long distances into Xqj bound states. 

Fig. 2. Color-octet diagram which mediates qq — > -i/'g production at 0{a1). 

Fig. 3. Lowest order contribution to the g* QQ^S^^] subamplitude in eqn. 2.5. 

Fig. 4. Color-octet diagrams which mediate gg — > "i/^Q production at 0{al). The sum of 
these graphs vanishes when both incoming gluons are on-shell. 

Fig. 5. Color-octet diagrams which mediate (a) qq i^gg, (b) gq if^gq and (c) gg — > 
ifjqg scattering at 0{al). The shaded circles appearing in these graphs represent 
the qq — > i/jq and gg i/jq amplitudes pictured in fig. ^ and fig. ^. 

Fig. 6. Theoretical transverse momentum differential cross section for T(IS') production 
at the Tevatron in the rapidity interval \y\ < 0.4 compared against preliminary 
CDF data. The dashed curve depicts the color-singlet contribution which includes 
direct T(15') production as well as radiative feeddown from Xbj{^P) and XbjC^P) 
states. The dot-dashed curve illustrates the color-octet contribution to T(IS') 
production. The solid curve equals the sum of the color-singlet and color-octet 
contributions and represents the total theoretical differential cross section. The 
T(15') cross section predicted by the high energy gluon fragmentation approxi- 
mation is shown by the dotted curve. 

Fig. 7. Theoretical transverse momentum differential cross section for T{2S) production 
at the Tevatron in the rapidity interval \y\ < 0.4 compared against preliminary 
CDF data. The curves in this figure are labeled the same as those in fig. ^ 
The dashed color-singlet cross section includes T(25') production and radiative 
feeddown from XbjC^P)- 

Fig. 8. Theoretical transverse momentum differential cross section for T(35') produc- 
tion at the Tevatron in the rapidity interval \y\ < 0.4 compared against pre- 
liminary CDF data. The curves in this figure are labeled the same as those in 
fig. ^. The dashed color-singlet cross section includes T(35') production and ra- 
diative feeddown from Xbj(3P). The measured n = 2 radiative branching ratios 
Br(xbj(2P) T{2S) + 7) have been used as estimates for their unknown n = 3 
counterparts Br(xbj(3P) — > T(3S') +7). 
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Fig. 9. Theoretical transverse momentum differential cross section for J/ij) production 
at the Tevatron in the pseudorapidity interval < 0.6 compared against prelim- 
inary CDF data. The dashed curve depicts the color-singlet contribution which 
includes direct J/V' production and radiative feeddown from Xcj states. The 
dot-dashed curve illustrates the color-octet contribution to J/?/' production. The 
solid curve equals the sum of the color-singlet and color-octet contributions and 
represents the total theoretical differential cross section. The J/ip cross section 
predicted by the high energy gluon fragmentation approximation is shown by the 
dotted curve. 

Fig. 10. Theoretical transverse momentum differential cross section for ip' production at 
the Tevatron in the pseudorapidity interval \ri\ < 0.6 compared against prelimi- 
nary CDF data. The curves in this figure are labeled the same as those in fig. ^ 
The dashed color-singlet cross section includes only direct i/j' production. 

Fig. 11. Ratio -Rtrans of transversely aligned and total QQ[^S[^^] differential cross sections 
plotted as a function of p±. The solid and dot-dashed curves illustrate -Rtrans for 
66['^5'^^-'] and cc[^5'^^-'] production respectively. 
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